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Abstract 

A magnetic system with a phase transition at temperature T c may exhibit double resonance peaks 
under a periodic external magnetic field because the time scale matches the external frequency at 
two different temperatures, one above T c and the other below T c . We study the double resonance 
phenomena for the mean-field g-state clock model based on the heat-bath-typed master equation. 
We find double peaks as observed in the kinetic Ising case (q = 2) for all q > 4, but for the three- 
state clock model (q = 3), the existence of double peaks is possible only above a certain external 
frequency since it undergoes a discontinuous phase transition. 
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I. INTRODUCTION 

From extensive investigations on stochastic resonance 3; it is now widely accepted that 
noise can play a constructive role. A schematic example of a particle in a double-well 
potential illustrates that the particle can move in a synchronized way with a weak external 
periodic force, when its average waiting time inside a well, determined by the noise strength, 
is comparable to the half of the period of the external forcing [1]. This is what is generally 
called a time-scale matching condition in studies of stochastic resonance. A particularly 
interesting model system for stochastic resonance is the kinetic Ising model since it can 
be regarded as coupled two-state oscillators with many degrees of freedom under thermal 



fluctuations [2J, y|. The probability to flip a spin j from Sj to — Sj is given by the Glauber 



aa. 

dynamics [4| as 

w(Sj -Sj) = - 1 - ^{.^(c^ v ; S A ) E(S U ...,-S s ,..., S N )]Y (1) 

where /3 is the inverse temperature defined by = 1/T with the Boltzmann constant ks = 1, 
and N is the system size. The energy E is a function of the spin configuration (Si, . . . , Sn), 
given by 

E = -jJ2 S i S i- h J2 Si ' 

(ij) i 

where J is a coupling strength, ^(y) runs over the nearest neighbors, and h is an external 
magnetic field. If we ignore spin correlations and assume that each spin experiences the mean 



field of the system, we arrive at the mean-field kinetic Ising model, w 



rich has served as an 
This model has 
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ideal starting point to study exact results on stochastic resonance 
been analytically shown that there can be two temperatures where the time-scale matching 
condition is met, one above the critical temperature T c of the system and the other below 
T c . The reason for this double stochastic resonance is that the intrinsic time scale of the 
system diverges in both the sides, whether T approaches T c from above or from below, so 
that the matching with the external frequency can occur on either side. Recently, the similar 
mechanism of time-scale matching in quantum kinetic Ising model is shown to be responsible 



for the double resonance peaks in quantum stochastic resonance Jfj]. 

A natural extension of the Ising model is the g-state clock model, where each spin 9 has 
an angle among discrete values 27m/ q where n = 0, . . . , q — 1. The spin at the jth site should 
now take a vector form as Sj = (cos 9j, sin 9j) = (cos(27m.,/g), sin(27mj/g)), and the energy 



function is accordingly rewritten as 

E = -J^2s i -S j -h'^2s i , 

(«) * 

where h = (h x , h y ) also takes a vector form. The total magnetization is M = N' 1 Si, and 
its magnitude M = \M\ will be used as an order parameter. The Ising model corresponds 
to the case of q = 2, and the XY model can be studied by taking the limit of q — > oo. If we 
wish to construct a kinetic dynamics for this g-state clock model, the probability to update 
spins such as Eq. (pQ) for the Ising case is readily obtainable by considering the heat-bath 
algorithm Q]. 

In this work, we check the double stochastic resonance in the g-state clock model within 
the linear response theory. When the external field h is parallel to the magnetization vector 
M, we find qualitatively the same double-resonance feature for all q > 3. When the field is 
perpendicular to M, however, the response is more complicated, and one of the resonance 
peaks will disappear as we approach the XY- model limit by taking q — > oo. We will pay 
particular attention to the case of q = 3, because the system undergoes a discontinuous 
phase transition unlike the other values of q > 2 {§(]. This work is organized as follows. 
In Sec. [Til we derive the equation of motion in terms of M for the general g-state clock 
model. From this equation of motion, we discuss the static phase transitions in Sec. IIHI 
Then, Sec. |IV] examines responses of the system when perturbed by a small amount from 
the static equilibrium. In Sec. |Vj we will see how the system under thermal fluctuations 
responds in a resonant way when the perturbation is given as a periodic magnetic field. 
Then, we conclude this work in Sec. I VI I 



II. EQUATION OF MOTION 

The master equation describing the probability distribution function P(0; t) for the spin 
configuration 6 at time t can be written as 

N N 

P(f):t) = 

dt 



±P(0; *) = "££ w 3 {6 3 -> tf.)P(9; t) + £ £ w,{9> -> 9,)P(0'; t), (2) 



where 6' = {9i, 9 2 , ■ • ■ , 9j, • ■ ■ , 9 N } differs from = 6 2 , ■ ■ ■ , 9j, ■ ■ ■ , 9 N } only at one site j. 
In the summations over 9L note that inclusion of the term for 9'a = 9j gives null contribution 
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in total, and thus has been included only for convenience. In the heat-bath algorithm, the 
transition rate is given by 

Wj{Bj ->• 0J) = ZJ 1 explPFj cos(^ - fa)), 

where 



with the inverse temperature (3. The local field is defined as 

V' = J ~ E ^ k + h ( 3 ) 

fceAj 

with magnitude and phase (f)j, where Aj is the set of nearest-neighboring sites of j and 
z is the coordination number (z = \Aj\ = const). We also denote the external local field 
hj as a complex number so that its real (imaginary) part yields the field in x (y) direction. 
The use of the heat-bath transition rate has a great benefit in calculation since it does not 
depend on the initial state, i.e., Wj(6j — > 6 1 -) = Wj{6'-), which enables us to write the master 
equation (j2J) as 

, N N 

j=l j=l $'. 



For an arbitrary single-spin function, denoted by f(9i), the following equation is derived 
from the master equation 

d , / E g ,exp^cos(g t -^)]/(ft) \ 

j7 /ft =" W + (-^ 7^ Tyi ) ) (5) 

at \ 2^ gi exp\pFicoa{9i-4>i)\ / 

as explained in Appendix |A] When q — 2, it recovers the kinetic Ising case in Ref. Q]. 
For the globally-coupled system with no external field, Eq. flH]) simply corresponds to the 
complex magnetization (we henceforth set J = 1) 

Me** = — el ° k ( 6 ) 

k=l 

and Eq. (j5J) turns out to be 

d (m)) _ (f(g)) , / E^exp^Mcos^-0)]/^) 



TABLE I. A qm = En=o cosm ( 27m /<7) for m = 0, 1, ■ • ■ ,4. We also list the critical value of the 
inverse temperature j3 c . For the g-state globally-coupled clock model /3 C = 2 for all q > 4. See text 
for details. 

q = 2 q=3 q=A q>5 

A q o q 
A ql 
A q2 2 q/2 

A q3 3/4 

A qA 2 9/8 4 3g/8 

(3 C 1 8 In 2/3 ps 1.848 39 2 

We further use /(0,) = e^"^ to get (/(0j)) = (1/JV) (E^"^) = ( M ) and 

d(M) / E 9; e^-^exp[/3Mcos(^-0)] \ 
* ^ ' \ E 0i exp[/3Mcos(^-0)] /• U 

The Hamiltonian of the g-state clock model without external field is invariant both under 
the uniform rotation, i.e., d\ — > Q\ + (ft, and under the reflection, i.e., Q\ — >■ —0;, which leads 
to 

d(M) / E ei cos^exp(/jMcos^) \ 
dt \ E ft exp(/3Mcos0,) /' U 

For the globally-coupled system in thermodynamic limit, the mean-field approximation be- 
comes exact and we can drop the average symbols to get 

dM = Eg, cosg; exp(/3Mcos0t) 

dt E e; exp(/3Mcos^) ' 1 j 

III. EQUILIBRIUM PHASE TRANSITION 

We are going to apply the above result in Sec. [II] to the g-state clock model, where 

Q\ = lunjq with n — 0, . . . , g — 1 in Eq. ffTO]) . In this section, we restrict ourselves to a static 

situation (dM/dt = 0) in the absence of the magnetic field (h — 0). In such a static case, 

Eq. (fTUj) is interpreted as = —dTjdM with the free energy J 7 , and assumes the following 

form of a self-consistent equation for M: 

M _ En=o cos(27rn/g) exp[/3M cos(27rn/g)] _ g In g 
En= 1 oexp[/3Mcos(27rn/g)] " d(/3M) ' 
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where 



q 1 00 (BM) m 

Z = J2 exp[/3M cos(27m/g)] = ]T ^J-A qm , (12) 

n=0 m=0 



with A qm = Xln=o cos m (2irn/q). Equation ( TTT|) is then expanded in the power of (3M as 

= Agi + ;.U.l, /2 + ^ 2 MM, 3 + |/3 3 M% 4 + Q(M 4 ) 
A q0 + (3MA q i + ^ 2 MM, 2 + \^M*A q3 + ^/3 4 M 4 + 0(M 5 ) 

= Aqi + A q pA q2 - A 2 ql + 2A 3 X - 3A q0 A ql A q2 + A 2 q0 A q3 
A q o A q0 2A 3 

6A* ql - 12^V4^y + AA 2 q0 A ql A q3 - gVV - 3^ 2 ] + 0(/?4m4) _ 

In Table [H we list values of A qm for m < 4. Since A g0 = g and A q \ — for g > 1, the above 
expansion is further simplified to 



.V = ^lm M ) + ^(/3M) 2 + qAqA 9 3A ' 2 (/3M) 3 + 0(/3 4 M 4 ). (13) 
g 2g eg"' 



The second term containing M 2 is particularly interesting, since it corresponds to the cubic 
term in J 7 , and thus is responsible for discontinuity of a phase transition [9j. Its coefficient 
in this case, A q3 , vanishes for every integer q > 1 except q = 3 (see Table H]). It agrees with 
our expectation since every mean-field g-state clock model undergoes a continuous transition 
except q = 3, which can be transformed to the mean- field three-state Potts model with a 



discontinuous transition 



A. q^3 



18 



When q 7^ 3, A q3 = (see Table H]) and Eq. (TIT]) does not have the cubic term: 



M a ^/3M + ^4 3 ^ (/3M)3 
g 6g 2 



From Table fl], v4 g2 = ^4 g4 = 2 for g = 2, and we find 

M « /3M- -/3 3 M 3 , 
3 

which yields M ~ (/3 — /3 C )^ 2 with /3 C = 1. For g = 4, we have A q2 = 2 and A q4 = 4, which 
yields 

2^ 24^ 
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FIG. 1. Magnetization of the q-state clock model as a function of T. (a) q = 4, 5, 6, and 7, which 
exhibit qualitatively the same behavior as M ~ \/3 — /3 C | 1//2 with /3 C = 2. (b) Three-state clock 
model (q = 3), exhibiting a discontinuous transition. The bistable region is bounded by 8i = 2 
and /?2 ~ 1.830 43 as depicted by the arrows. The vertical dashed line indicates the transition 
point C 



| In 2 » 1.848 39. 



resulting in M ~ {8 - (3 C ) 1/2 with (3 C = 2. This agrees with an exact relationship between 
q = 2 and 4 



10|. For q > 4, we have A 92 = q/2 and ^4 9 4 = 3g/8 and we arrive at 
M 



-BM - —B 3 M 3 , 
2' 16 



irrespective of q, which means that we always find M ~ (8 - d c ) 1 / 2 with C = 2 [Fig. QJa)]. 
The scaling form of M with respect to the temperature also confirms the mean-field value 
1/2 of the magnetization critical exponent, for all values of q other than three. We finally 
remark that for the XY model (q — > oo), we replace the summations in Eq. ( II ip by integrals 

as 

' dde Wco S e\ h{BM) 



M 



d 



d(8M) 



In 



where I„. is the modified Bessel function Uj. It is known that 8 C 



h{8M) 



(14) 

2 in the limit of 



q — > oo [III, ll2j, which is in agreement with the above conclusion of 8 C = 2 for q > 4. These 
are verified by numerical calculations as shown in Fig. [TJa), which are obtained by directly 
solving Eq. ([11]) in a numerical way. 



B. q = 3 

In the three-state clock model, we observe a discontinuous phase transition due to the 
cubic term in the expansion of the self-consistent equation (fTTj) . However, in order to look 
into the bistable region in detail, we are not allowed to use the expansion in terms of 8M, 
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FIG. 2. Free-energy landscapes of the three-state clock model as a function of M. The free energy 
T = M(M + l)/2 - log(e 3/3M / 2 + 2) is obtained from the integration of Eq. 1(15]). (a) If 
/?>/?! = 2, there exists only one free-energy minimum at nonzero M. If /? is decreased further 
below Pi, T begins to have two minima, but the global minimum of T occurs at nonzero M for 
(3 > f3 c . (b) At (3 = j3 c = | In 2 ~ 1.848 39, the two free-energy minima have equal height, one at 
M = and the other at M = 1/2. As /3 is decreased below /3 C , the minimum at M = begins 
to locate lower than the other minimum at nonzero M, but the two minima coexist until P = P2 
is reached, (c) When j3 < ~ 1-830 43, the only minimum is found at M = 0. Note that the 
system is bistable for /?2 < f3 < (3\. The vertical dashed lines indicate free-energy minima, and the 
values of M at these minima are observed in Fig. QJb). 

since M cannot be assumed to be small in this case. We thus start from Eq. ( TTUj) with the 
equilibrium magnetization M which should satisfy 

e 3/3M/2 _ 1 

= -M+ eWM/2 + 2 - (15) 



Let us assume that there exists a certain (3, where the derivative of Eq. ( 1151) with respect 
to M vanishes. This tells us when bistability becomes possible. Denning y = e 3/3M//2 , these 
two conditions can be written as 

y + 2' 
/3 = 2^ 1 (y + 2) 2 /9, 
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which leads to 

1 + 2M 
V ~ 1-M' 

2 

^ = (1 + 2M)(1-M)' ^ 
Inserting these into the definition of y, we get a transcendental equation for M as follows: 

1 'nfV^V W 



1-M 1 + 2M V 1 - M 

Note that there is a trivial solution M\ = with fli — 2 from Eq. (fl6l) . One can also 
find a nontrivial solution of Eq. (fl7|) numerically as M 2 ~ 0.377 201, and the corresponding 
inverse temperature is 02 ~ 1.830 43. It means that the system is bistable between /3i and 
[Fig. GJb)]. The transition point (3 C can be determined by checking when the two free-energy 
minima have an equal height as in the Maxwell construction. Interpreting the left-hand side 
of Eq. (T15D as dT/dM, we find that EL = fin 2 » 1.848 39 (Fig. [2}, which is between ft 
and /?2 and in agreement with Refs. One can also readily check that the nonzero 

magnetization at (3 C is M = 1/2 by using Eq. (p~5|) . 



IV. RELAXATION TIME 

Returning back to the general g-state case, the evolution equation of the magnetization 
vector (M x , M y ) is given by putting f{6i) = cos Q\ and f(Q\) = sin 9i in Eq. §7§ since 
M x = (cos 9i) and M y = (sin #/). Within the mean-field scheme with no external field, we 
find 

dM. = _ M E g cosflexp[/3Mcos(fl-0)] 

dt x £ e exp[/3Mcos(fl -</>)] 

dM y = E e singexp[/3Mcos(g-0)] 

dt " y ^2 exp[/3Mcos(9-(f))} ' 

where (M x ,M y ) = (M cos 0, M sin 0) from Eq. ([6]) and runs over = 2irn/q with n = 
1, . . . , q — 1. We then add perturbation 6 = (8 X , 5 y ) around the equilibrium magnetization 
with the assumption /3\d\ <C 1. By expanding the exponential functions and leaving linear 
terms with respect to (35 x or (35 y , one can study linear responses of the system. Without loss 
of generality, we may set the initial magnetization along the x-axis, i.e., M* = (M*, M*) = 

9 



(M*,0), where M* satisfies Eq. ([II]) as 

= EflCosflexpQgM* cosfl) 
£ e exp(/3M*cos0) ' 
_ ^ sin 9 exp(/3M* cos 0) 
~~ ^ 

The linearized equations for perturbation S x (t) = M x {t) — M* and 8 y (t) = M y (t) — M* are 
derived in Appendix [B] as follows: 

^ = -{l + (3M* 2 -f3C)5 x , (18) 
^ = -(l-)9 + 0CH, (19) 

where 

E e cos 2 flexpQgM* cosfl) 

£ e exp(/3M*cosfl) ' 1 J 

In other words, we have two relaxation times rii and rj_ in the parallel and perpendicular 

directions with respect to the equilibrium magnetization vector, respectively, which are given 

as 

t- 1 = 1 + f3M* 2 - PC, (21) 
rl l = 1 - (3 + f3C = 2 + f3(M* 2 -1) -t^ 1 . (22) 

Very near to the critical point, we may assume that M* « 0, and we get C = ^ e cos 2 6/ ^ e 1 = 
A q2 /A q0 (see Tabled]). When q = 2, the y component of magnetization is not defined, and 
we have C = A22/A20 — 1 from Table H] yielding t^ 1 = 1 — j3. For q > 2, we instead 
have C — 1/2 from Tabled] and Ty" 1 = r^ 1 = 1 — /3/2. Consequently, the divergence of 
the relaxation time occurs when = 1 for q = 2, and when (3 = 2 for g > 2 (g 7^ 3), in 
complete agreement with (3 C = 1 for q = 2 and /3 C = 2 for g > 2 (except for g = 3 in which 
M* rs is not justified) found in Sec. IIII Al Also, the divergence of the relaxation time in 
the form of r ~ |/3 — /3 C | _1 , determines the dynamic critical exponent z = 2 in the scaling 
form r ~ £ 2 , since the correlation-length (£) exponent v = 1/2 in £ ~ |/3 — /3 C |~ ! '' for the 
mean-field universality class. 



A. 



As mentioned above, t± becomes meaningless for the Ising case (q = 2), and we get 
- 1 = l + (3M* 2 -(3, recovering the result in Refs. [Fig. IHJa)]. For different values of 
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FIG. 3. Inverse relaxation times for the g-state clock models, (a) q = 2, (b) q = 4, (c) q = 5, and 
(d) the XY limit — > oo) where r7 — > above /3 C = 2. 



g, we plot Eqs. (12T|) and ( 122|) in Fig. [3] by solving Eq. (TTTj) . For q = 4, we find ry = in 
the entire region of T [Fig. [3]^b)] since 

2 



c 



e PM* _|_ e -/3M* 
e /3M* + e -/3M* + 2 



1 + 



;^ M * + 1 



1 + M* 



For q > 4, we instead observe t_i_ > rn at (3 > C [Fig. E(c)] and the difference becomes 
pronounced as q increases. In the limiting case of q — > oo, t± diverges at (3 > (3 C , which 
reflects the U(l) symmetry of the system. This shows the following identity of the modified 
Bessel function for (3 > 2, 



2 



1 + 



WM* 



1 

'-J 



WM*)\ 

where M* satisfies Eq. f lT4|) . This observation is also related to the idea in Ref. 14| that 
fluctuations in the angular direction can distinguish the discrete symmetry in the clock 
model from the continuous symmetry of the XY model. 



B. q = 3 



As explained above, this system has a bistable region between Pi = 2 and p2 ~ 1-830 43. 
Even though the relaxation times can be expressed in the same way, one should note that 

11 
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FIG. 4. Three-state clock model, (a) Inverse relaxation time t7 in the direction of magnetization 
and (b) rj 1 in the perpendicular direction. The solid lines represent values at the lowest free 
energy, while the dotted lines mean what one can observe at metastable states. The vertical 
dashed lines indicate the transition point j3 c = | In 2 1.848 39. The horizontal dashed line shows 
1- |ln2 » 0.537 902. 

they are defined in the linear-response regime. In other words, even if the system becomes 
bistable, the relaxation here means returning back to the original state and not jumping to 
the other state. If M* = 0, we have seen that r^ 1 = r^ 1 = 1-/3/2. So both of them vanish 
if /3 — > fli from below with keeping M* = 0. If /3 approaches (3 2 from above with keeping 
M* ytz o, on the other hand, we find 

-I fa , a+\ 1 ^ 2 . fl U 2 n 

r H (P = 1 - y y~ + ^ 2 2 = ' 

by using Eq. (|T6|) . Here M 2 means the nontrivial solution of Eq. (IT7|) . These results are 
plotted in Fig. H](a). One can clearly see why those metastable branches cannot be sustained 
beyond f5\ and fa, respectively. The other inverse relaxation time rj 1 is plotted in Fig.@Jb). 
Note the jump at & from rj 1 ^ -»■ = 1 - /3 c /2 = 1 - | In 2 « 0.075 803 8 to tJ 1 ^ -»- 
= 1- | In 2 » 0.537 902. The latter value is obtained by using Eq. ((22} with M* = 1/2. 
This behavior is in accordance with the general tendency that t± increases in the ordered 
phases as q becomes larger, but manifests itself in a discontinuous way. 



V. RESONANCE 



If there exists a uniform external field h = (h x , h y ), the evolution equations are generalized 

to 

dM x E 9 cosgexp[/3Fcos(g-0)] 
dt x £ e exp[/3Fcos(#-0)] 1 ' 
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dMy E e sinflexp[/3Fcos(fl-0)] 
dt v £,exp[/3Fcos(#-0)] ' 1 ' 

where (Fcos0, Fsin</>) = (M x +h x , M y + h y ) [see Eq. (j3J) for globally-coupled (2; = iV) model 
with a uniform external field]. The magnetization is again decomposed into (M x ,M y ) = 
(M* + 8 X , M* + 8y), where we set again M* = M* and M* = as before. Let us furthermore 
assume that (3\h\ <C 1 and expand the equations up to the linear order. The derivation is 
almost the same as what we did for Eqs. (fl8l) and (fT9l) . except that we have to replace S x by 
5 X + h x as well as S y by 5 y + h y in expanding the second terms of Eqs. ff23|) and (JM]) . This 
results in 

^ = -5, - ft (M* 2 - C) (5. + M 

= — T-jj -1 ( Ja; + h x ) + 

= -r^8 x + (1 - r- 1 )^., 
^ = -^ + /3(l-C) + 



— T A 1 {5y + Hy) + k 



y 



where t^ 1 and rj 1 are given by Eqs. (}2~Tj) and (1221) . Since these equations have basically 
the same form, we may drop the subscript || or ± in finding the formal solution. If we set 
h = h cos ut, the resulting equation is the following: 

dS 

— = -r _1 <5 + (1 -r^hocosut. (25) 

(JjL 

Its solution is obtained by assuming 

5(t) = 5o cos(wt — a), 

which results in phase shift a = arctan(wr) and amplitude 

(r-l)h 



So 



2^2 



We can interpret this solution as follows. When r ~ 0(1), the left-hand side of Eq. ( I25i) 
is negligible compared to the first term of the right-hand side since 00 <C 1 by assumption. 
Then, 5 becomes directly proportional to ho cos cut without any phase shift, even though the 
amplitude will be small. On the other hand, if r 3> 1, the first term on the right-hand side 
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of Eq. ( )25|) becomes negligible so that dS/dt is proportional to hocosut, which leads to a 
phase shift of 7r/2 in 5(t) compared to the magnetic field. 
The ac susceptibility components are defined as 



1 



2tt 



X = —— / d(ut) Mcoswt, 

TT^O Jo 

i r 2n 

X = —— / d{ut) Msmut. 

7™0 Jo 

I 1 

The former one x' is closely related to the occupancy ratio used in Ref. [3| to measure how 
many spins are aligned in the direction of the external field. By inserting S(t) here, we 
obtain 

*' = TT^' (26) 

X" = = (27) 

where the static value M* integrates out to zero when multiplied by the sinusoidal functions. 
When M* = and r" 1 = 1 - fi/2, the maximum of Eq. Q2SJ) is found at /3 = 2w 2 + 2 - 
2uj\/u} 2 + 1, which approaches unity as u — > oo. In general, the extremum condition of x' 
with respect to (3 is equivalent to 2o;x" = 1 since 

d X ' dr/dp t-1 2 dr / dr\ 1 - 2uj X " 

~dj ~ 1 + uV 2 ~ (l + w 2 r 2 ) 2 ^ T rf/3 ~ V W 1 + oo 2 t 2 ' 

This condition yields a solution w = (r 2 — 2r) -1 / 2 , which can be expanded as r^ 1 + r~ 2 + 
|r~ 3 +- • •. Therefore, the optimal frequency u for resonance coincides with r^ 1 to the leading 
order when r 3> 1. Consequently, the stochastic resonance occurs when the extrinsic time 
scale 1/u; matches with the intrinsic one r 



A. q^3 

First, we consider the field in the x direction, parallel to the magnetization. We thus use 
t\\ in place of r in Eqs. (|26l) and ( 1271) . If the system undergoes a continuous phase transition 
with g^3, one can clearly see two peaks in x' above and below /3 C . We show the case of 
q = 6 in Fig. GD^a), noting that qualitatively the same behavior is observed for other q values. 
Now let us apply the field in the y direction for q ^ 2 (note that the field in y direction 
for q = 2 is meaningless). When (3 < C , the disordered phase of the system is isotropic so 
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FIG. 5. ac susceptibility components for q = 6 given by Eqs. (|26p and (|27p with oj = 27r x 10 _1 . 
(a) The field can be exerted in the x direction or (b) in the y direction. The next two panels show 
the same plots for q — ¥ oo, with (c) the field in the x direction and (d) in the y direction. The 
horizontal dashed lines indicate x" = (2a) _1 to locate the maxima in x' ■ 



we observe the same response as above, although the field direction has changed. However, 
when (3 > (3 C , the peak is suppressed to a higher (3 [Fig. [5](b)]. As q — > oo, x' eventually 
vanishes and x" becomes constant at > (3 C , which signals the U(l) symmetry [Fig. 0(d)] . 



B. q = 3 

When the field is parallel to the magnetization, one can infer from Fig. IU(a) the possibility 
of two resonance peaks only above a certain frequency since one can find two temperatures 
where the external frequency matches the relaxation time scale only if r _1 > 1 — c /2 ~ 
0.075 803 8. One cannot find proper relaxation time scales to match the external driving 
frequency if the frequency is too low. It is readily confirmed in Figs. [3(a) to [6](d). More 
precisely, the threshold of u for the stable double resonance peaks can be found by requiring 
X" = (2a;)" 1 to be met exactly at (3 = (3 C = § In 2. Since r" 1 ^) = 1 - (3 c /2, we get the 
following quadratic equation, 

WT||(T[|-1) _ 1 



l + w 2 r,2 2a; 5 



(28) 



which yields u = (3 - 4 In 2) In 2 - 9 w 0.0822986. 
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FIG. 6. Three-state clock model under the field in the x direction, (a)-(b) ac susceptibilities at 
to = 2ir x 10 _1 and (c)-(d) the same quantities but at oj = 2ir x 10~ 2 . The solid lines represent values 
at the lowest free energy, while the dotted lines mean what one can observe at metastable states. 
The vertical dashed lines indicate the transition point /3 C = §hi2 ~ 1.848 39. The horizontal 
dashed lines indicate (2a;) . 

When the field is perpendicular to the magnetization, one can guess that the external 
frequency u should be again large enough to find two matching temperatures. In Fig. H^b), 
for example, the minimum tJ 1 above which there stably exist two matching temperatures 
is found to be rj 1 = 1 — |ln2 « 0.537 902. It is true that one finds only one maximum 
in x' when u is low [Figs. [7Ja) and (b)]. However, our calculation shows that the double 
resonance is anyway impossible if we take only truly stable states into account. The left 
maximum in x' can De located only on a metastable branch even for a very large value 
of uj [Figs. [7(c)] . If we repeat the same calculation as Eq. ()28l) to find the threshold of uj 
with tJ 1 = 1 — I In 2, we indeed find that the equation does not possess any real solution, 
confirming this impossibility. 



VI. SUMMARY 



In summary, we have studied stochastic resonance with the mean-field kinetic version 
of the g-state clock model. The response under a periodic external field now depends on 
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FIG. 7. Three-state clock model under the field in the y direction, (a)-(b) ac susceptibilities at 
lo = 2tt x 10 _1 and (c)-(d) the same quantities but at w = 2ir. Inset: a zoomed view of the left 
branch, showing a maximum around /3 ~ 1.837. The solid lines represent values at the lowest 
free energy, while the dotted lines mean what one can observe at metastable states. The vertical 
dashed lines indicate the transition point (3 C = |ln2 ~ 1.848 39. The horizontal dashed lines 
indicate {2uj)~ 1 . 



the direction of the field relative to the magnetization vector. When they are parallel, the 
double stochastic resonance is observed for every q > 3 qualitatively in the same way as in the 
kinetic Ising case (q — 2) jjj. When the field is perpendicular to the magnetization vector, 
on the other hand, the resonance peak is suppressed to a lower temperature and eventually 
vanishes as q — > oo since the U(l) symmetry sets in. For q = 3, the discontinuous transition 
should be also taken into account, and we have concluded that the double resonance peaks 
are observable in a truly stable manner only when the external driving frequency is high 
enough and the field direction is parallel to that of the magnetization vector. 
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Appendix A: Derivation of Eq. ([5]) 



We multiply each side of Eq. @ by an arbitrary function of spin I, denoted by f(9i), and 
carry out summation over all the possible configurations: 

, N N 

E j t p ^ = - EE + E E E w M -+ °i) p w ( A1 ) 

e e j=i e i=i ^ 

The sum over j is decomposed into two parts: the sum over terms with j ^ I and the term 
for j — I. Let us consider the former first. If we write 9j = 9 3 ■ + A, we get 

N N 

"EE P ^ WW + E E E w iVj) p ( e i> ■ ■ ■ ,9j + A, - ■ ■ ,9 N ; t)f{9 l ) 
j+x e j^i a 

N N 

= - E E p (0- o/w + E E E w M - = o, 

e jjti e &i a 

where £\, iw,-(0i-) = 1 has again been used. With (/(0,)) = Ee P(0; t)/(0,) and 0,' = 0, + A 
again, Eq. flAlj) is now reduced to 

j (f(9i)) = ~ <M)> + E E w i(0i)P(0n ■ ■ ■ , 9 l + A, • • • , 9 N] t)f{9 l ) 

e a 

= - (M)) + E E w ^ - A ) p ( - - A ) 

A 

= -(/(0z)> + (E^- A )^- A ) 

= " (/(ft)) + 

= - (M)) + 



A 

E A exp[/3F,cos(ft A)]/(g ; - A) 

Ee, exp[/3Fi cos(0j - </>;)] 
E e , exp^cos^-^)]/^) 



Ee, exp[/3F z cos(6» ; - 0,)] 
Appendix B: Derivation of Eqs. (I18|) and fjl9j> 

The perturbation in M x is expanded up to the first order as follows: 

d5 J2 CQS Q e PM, cos 6 e PM y sin 9 

- -Ml - 5 X + ^ 6 



x x ^2 eP M x cos O e fiM y sine 

ECOS/9 e PM* coaOgPSxCosOgPSyamt 
6 



-M* -5 X + 
-M* -5 X + 



e /3Af * cos eg/S^a; cos 6» e /3(5 y sin 

E e cos 9 e? M * cos8 (l + (35 x cos 9 + (35 y sin i 
E e e^* cos9 (l + Wx cos + ^ sin 9) 
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_ M * Eecosfl e P M * cose + J2 9 cos6 e m * cose {f35 x cos 9 + (35 y sin 9) 

x+ Eg e^ M * cose + e^ M * cose (/?4 cos 9 + ^ sin 9) 

« -M* - 4+ < J^costf e^ M * cose + ^cos# e^ M * cose (/?4 cos + sin 0) 

. e e J I e 

= ~ 5 *~ \j2 cose e^ cos ^e^ cose (^cos^ + ^sin^)} x <j 



D /3M* cos6> 



2 

COS0 

-2 



+ J^e /3M * cos9 ^cos^e /3M * cose (/34cos^ + /3^sin^)l x i J^V M * 

I J I 

-5 -/3M* 2 5 Eg cos exp(/3M* cos 0) (ffi a cos 9 + (55 y sin g) 
x ' E e exp(/3M*cos0) 

Likewise, we obtain 

dS v f E0 sm ex P(£ M *cos0)(/^cos0 + /^sinfl) 
"dt" ~ ~ y + ^ exp(/3M* costf) ' 

One can furthermore show that 

^2 cos 9 sin 61 exp(/3M* cos 6>) = 0, 
e 

by which we can rewrite the above equations as 



cos 



where 



dt 

d6y - = - (I - (3 + (3C) 5 y , 



c _ Eg c °s 2 exp(/3M* cos 0) 
= J2 d exp(/3M* cos 9) 
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